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7^ I Abstract 

I . The shifted - / expansion technique (SLET) is extended to solve 

"^ I for Dirac particle trapped in spherically symmetric scalar and/or 4- 

^ ' vector potentials. A parameter A = 0, 1 is introduced in such a way 

that one can obtain the Klein-Gordon (KG) bound states from Dirac 

bound states. The 4- vector Coulomb, the scalar linear, and the equally 

k> I mixed scalar and 4-vector power-law potentials are used in KG and 

j_j ■ Dirac equations. Exact numerical results are obtained for the 4-vector 



Coulomb potential in both KG and Dirac equations. Highly accurate 
and fast converging results are obtained for the scalar linear and the 
equally mixed scalar and 4-vector power-law potentials. 



I Introduction 

In a previous paper we have developed the shifted-/ expansion technique 
(SLET) to solve 3D and 2D Schrodinger equations [1]. Therein, we have 
suggested SLET as a reformation to the shifted large-N expansion technique 
(SLNT) [2,3]. The importance of the Klein-Gordon (KG) and Dirac wave 
equations with scalar and/or 4-vector potentials arises in many fields of 
physics [4-24]. In the framework of the mentioned equations, a scalar or 
an equally mixed scalar and 4-vector potentials have considerable interest in 
the study of quarkonium systems [6,8-10,12,16-19]. The mixture of 4-vector 
and scalar potentials has been studied by Long and Robson [20]. The 4- 
vector potentials have great utility in atomic, nuclear, and plasma physics 
[21-24]. Therefore, many attempts have been made to develop approximation 
techniques to treat relativistic particles in the KG and Dirac wave equations 
[4-24]. 

To the best of our knowledge, the attempts that have been made to 
develop a relativistic 1/A^ expansion technique include, Neito [14] who has 
extended the unshifted 1/A^ expansion [25] formalism to treat a vr-mesonic 
atom in the KG equation. Miramontes and Pajares [15] who have studied the 
same technique in the framework of the KG and Dirac equations. Iterative 
procedures concerning the 1/N expansion technique were introduced to deal 
with KG and Dirac particles [11]. Stepanove and Tutik [10] who have used 
the /i-expansion procedure in a similar manner to the 1/A^ expansion proce- 
dure. Such procedures are straight forward though rather tedious and suffer 
from convergence problems in comparison to SLNT [4-8]. In addition, highly 
excited states pose problems for the above mentioned techniques. Recently, 
Roychoudhury and Varshni have developed SLNT to treat scalar [12] and 
4-vector potentials [13] in the Dirac equation. Panja et al [4] have applied 
SLNT to Dirac particle by considering the KG equation in which a spin-orbit 
interaction is introduced analogous to Pauli theory ( this has been proposed 
by Papp [8]). Mustafa and Sever [6-8] have used a different approach to 
SLNT for the KG and Dirac equations. 

Encouraged by the success of SLET in the Schrodinger equation and the 
importance of KG and Dirac equations with scalar and/or 4-vector potentials, 
we feel tempted to extend SLET to solve for the bound states of the above 
equations. We shall do so by considering the Dirac equation in which we 
introduce a parameter A = 0, 1 in the spin-dependent term. For A = 0, 



Dirac equation implies KG equation with scalar and/or 4-vector potentials. 
In this case, the total angular momentum j is given through the relation 
I = j ± 1/2, and the radial wave function of KG particle is the radial large- 
component of the Dirac spinor. When A = 1, the following are the bound 
states under consideration; (i) the Dirac bound states in scalar and/or 4- 
vector potentials, and (ii) the Dirac bound states in an equally mixed scalar 
and 4-vector potentials. The bound states in (ii) are identical to KG bound 
states in the same structure potentials. 

For the sake of numerical illustration on the accuracy of our SLET we 
discuss the results for the following particles; (i) Dirac and KG particles 
in 4-vector Coulomb potentials, (ii) Dirac and KG particles in scalar linear 
potentials, and (iii) Dirac particle in an equally mixed scalar and 4-vector 
power-low potential. Elsewhere, we shall investigate Dirac and KG particles 
in a mixture of 4-vector and scalar potentials. 

In Sec II we shall extend SLET to solve for Dirac and KG bound states. 
Analytical expressions are given in such away that allows the reader to use 
them without proceeding into their derivations. In Sec III we shall show that 
the analytical expressions of Sec II lead to highly accurate and fast converging 
numerical results when applied to the 4-vector Coulomb, the linear scalar, 
and the equally mixed scalar and 4-vector power-law potentials in both Dirac 
and KG equations. We conclude and remark in Sec IV. 

II SLET for Dirac and KG bound states. 

With a scalar S{r) and a 4-vector V{r) Lorentz interactions, simultaneously 
present, the radial coupled Dirac equations can be written [16,26] {h = c = 1) 
as 

eiG(r) + ^ - -Fir) = 0, (1) 

ar r 

6F(r) _ ^ _ tcir) = 0. (2) 

ar r 

Where ^i = E — m — V{r-) — S{r), ^2 = E + m — y{r), y{r) = V{r) — S{r), 
and E is the relativistic energy. G{r) and F{r) are the large and small radial 
components of the Dirac spinor [19], respectively. Provided V{r) and S{r) 
are spherically symmetric, k = — (/ -|- 1) for the total angular momentum 



j = I + 1/2, k = I for j = I — 1/2, and / is angular momentum quantum 
number. 

In terms of the large-component G{r) of the Dirac spinor, Eq.(l) reads 






r^ ^2 df ^ 



G(r) = 0, 



(3) 



where the prime denotes d/dr. 

For any k one can show that k{k + l) = /(/ + 1). Furthermore, we remove 
the first derivative by proposing the ansatz 



G{r) = R{r)exp{—p{r) /2) ; p (r) = y (r)/^2- 
Which in turn implies the Schrodinger-like radial Dirac equation 

-d^dr^ + l{l + l)/r2 + f/(r) - ^61 Rij') = 0, 



with 



U{r) 



1 
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\2y"{r) 
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-'^^'(^)+3p(rf] 



(4) 



(5) 



(6) 



It should be pointed out that the choice of the ansatz in Eq.(4) is mo- 
tivated by requiring the agreement between Eq.(5) and the outcome of the 
linear matrix eigenvalue problem, Eq.(l) and (2), after being transformed 
into a diagonal Sturm-Liouville eigenvalue system. For more details on this 
system the reader may refer to Barut [27]. Moreover, it is evident that Eq.(5) 
reduces to KG wave equation with scalar and/or 4- vector potentials provided 
U{r) = [6]. We therefore introduce a parameter A = 0, 1 in U{r) so that 
A = and A = 1 correspond to KG and Dirac wave equations, respectively. 
Also, we shall be interested in problems where the rest energy is large com- 
pared to the binding energy W = E — m. Eq.(6) thus becomes 



f/(r) 

where we have considered 
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2m 



W - y{r) 
2m 



+ 



(7) 
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and the term {W — y{r))/4m'^ is small enough to be neglected. 
If we shift the angular momentum / through the relation / 
define 



7(r) = -V{rY + m{rY + f/(r), 



where 



m{r) =171 + S{r) 
Eq.(5) becomes 

d^ [P + li2/3+l)+(3i(3 + l)] 



dr'' 



+ 



+ 7(r) + 2E\/(r) 



We shall now start the systematic 1// expansion by defining 
p 
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E, + Ei/l + E2/I' + E^/l' + 



/ + j3 and 
(9) 

(10) 

R{r) = E^R{r). 

(11) 

(12) 
(13) 
(14) 



where x = l^^'^{r — To)/ To, and Q is to be set equal to t^ at the end of the 
calculations. With Eqs.(12)-(14) into Eq.(ll) one gets 
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Eq(15) is a Schrodinger-like equation for the one-dimensional anharmonic 
oscillator and has been discussed in detail by Imbo et al [2]. We therefore 
quote only the resulting eigenvalue of Ref.[2] and write 
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where ai and 0:2 are given in the appendix of Ref. [1], and the corresponding 
relativistic e s and 5 s are given in the appendix of this paper. If we now 
compare Eq.(17) with (16), we obtain 
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To is chosen to be the minimum of Eo, i. e.; 



dEo/dr-o = and d'^Eo/drl > 0. 

Hence, Tq is obtained through the relation 

2(/ - pf = b{ro) + ^b{roy-4c{ro), 



(23) 



Where 



b{ro) = rl 2Viro)V\ro) + t'I^o) + ToV' (rof 



c{ro) = -f 



(24) 

(25) 
(26) 



7 {rof + W{ro)V {roh (^o) - ^l{ro)V {ro)\ 
The shifting parameter f3 is determined by requiring £"1 = [1,2]. Therefore, 

P = -[l + {nr + l/2)w]/2, (27) 



where 
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III Applications, numerical results, and dis- 
cussion. 

In this section we shall consider the 4-vector Coulomb, the scalar linear, 
and the equally-mixed scalar and 4-vector power-law potentials in KG and 
Dirac equations. To obtain the KG and Dirac bound states we set A = 
and A = 1, respectively. Tq is obtained from Eq.(24) through (25)-(29), and 
the energy eigenvalues are calculated from Eq.(22). We shall examine the 
accuracy and the convergence of SLET by comparing its outcomes with the 
exact and numerical results of other methods. 



111.1 The 4- vector Coulomb potential in KG and Dirac 
equations [V{r) = —a/r, S{r) = 0] 

A pionic atom in a Coulomb-type 4- vector potential, V{r) = —a/r and 
S{r) = 0, obeys the KG equation. Where a is the fine structure constant 
and a = 1/137.03602. The pion mass is m-irC^ = 139.577 MeV and spin-0. 
An electron trapped in the same type of potential obeys the Dirac equation. 
The rest mass of the electron is meC^ = 0.5110041 MeV. However, for both 
particles exact analytical results for binding energies exist [26] and hence 
exact numerical results can be reproduced and compared with those of SLET. 
In Tables 1 and 2 we show the results of SLET compared with the exact 
results for KG and Dirac particles, respectively. The results are shown in 
such a way that one can easily judge the accuracy and the convergence. Our 
results appear in excellent agreement with the exact ones. The minimum 
and maximum percent accuracies found for both particles are 99.998% and 
100.000%, respectively. The convergence is noted to be very fast in the sense 
that the leading term Eg contributes from 99.996% to 100.000% of the exact 
binding energy. 

111. 2 The scalar linear potential in KG and Dirac 
equations [V{r) = 0, and S{r) = Ar] 

A pure scalar linear potential, V{r) = and S{r) = Ar, in KG and Dirac 
equations is precisely a quark confining linear potential. This potential has 
been studied by Gunion and Li [17] in KG and Dirac equations and hence 
accurate numerical results exist for comparison with SLET results. However, 
Mustafa et al [6] and Roy et al [12] have used SLNT to calculate part of the 
KG and Dirac mass spectra, respectively, of some quark- antiquark systems. 
In Tables 3-5 we compare SLET results with those of SLNT [6,12] and 
those of Gunion and Li [17]. Comparing our results, in Table 3, with the nu- 
merically predicted ones [17] we have scored a minimum accuracy of 99.96% 
and a maximum accuracy of 100.00%. In Table 4 the minimum accuracy is 
noted to be 99.61% and the maximum is 99.76%. In Table 5 the accuracy 
ranges from 99.05% to 100.00%. Furthermore, the energy eigenvalues given 
by Eq.(14) have been noted to be fast converging in the sense that the leading 
term Eo contributes more than 99% of the total energy. While our results in 
Table 3 are in exact agreement with those of SLNT [6] , they are only in qual- 



itative agreement with those of SLNT [12] in Tables 4 and 5. The term that 
appears as 3y {rY/Sm^ in Eq. (7) of this text has appeared as y (r)^/8m^ 
in Eq.(lOa) of Ref.[12]. Their method [12] was restricted only to the scalar 
linear potential S{r) = Ar in the Dirac equation where the effective poten- 
tial can easily be treated by the usual nonrelativistic SLNT. In addition the 
authors of Ref.[12] were interested in terms only up to 0{A/m^), a thing we 
have not assumed here. However, the term 3y [rY/Sm'^ has been confirmed 
to exist by Barut [27]. Moreover, it is noteworthy to mention that for some 
values of n^ and / the standard numerical method [17] could not predict the 
effect of the spin-orbit coupling term L.S. SLET was sensitive enough to 
predict this effect. To see this one may compare the results of Table 3 with 
those of Table 4. 

III. 3 The equally mixed scalar and 4-vector power-law 
potential [V{r) = S{r) = Ar" + K] 

The power-law potential of the form 

V{r) = Ar" + Vo (29) 

with z/ = 0.1 and A > is a simple non-QCD- based potential. An equally 
mixed scalar and 4-vector structure of this potential, i.e. V{r) = S{r), in 
Dirac equation is known to reproduce the data of \E' and T spectroscopies. 

In Tables 6,7,8 and 9 we show the spin-averaged results for the equally 
mixed scalar and 4-vector power-law potential, Eq.(29), with Vo = —2.028 
GeV, u = 0.1, A = a''+\ a=1.709, m^ = 1.6179 GeV, rub = 5.0114 GeV, 
rris = 0.325 GeV, and ma = 0.01 GeV. 

In Table 6 we compare the eigenvalues e„^/ corresponding to the confined 
bound states of quarks. Sn^i is calculated through the relation [19] 

r o / iu/(v+2) 

en^i = iE,-m,-2Vo)[iE, + m,) {2 A) ^^Z-] (30) 

In this table the accuracies of our results are in the range 99.90% to 100.00%. 
Whereas in Ref.[5], SLNT has scored a maximum error of 2.7% correspond- 
ing to a minimum accuracy of 97.3%. Therein, the authors [5] should have 
noticed that the term on the right hand side of Eq.(5), in Ref.[5], vanishes for 
the equally mixed scalar and 4-vector potential case. This makes it unneces- 
sary to seek the nonrelativistic limit they have sought. The consideration of 
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this limit in Ref . [5] has caused a spin- orbit couphng term to exist (see Eq. (9) 
of Ref. [5]) although it is well known that one of the interesting features of 
such mixture cancels any spin-orbit coupling effect [18]. In Tables 7, 8 and 
9 we investigate the accuracy of SLET for some heavy, bb, light, ss, and 
non-self-conjugate atom like cd mesons, respectively. Therein, SLET results 
are found comparable to the numerically predicted ones [19]. 

IV Conclusions and remarks 

In this paper we have extended the shifted-/ expansion technique (SLET) 
to solve for the eigenvalues of KG and Dirac equations with scalar and/or 
4-vector potentials. We have shown that Dirac equation reduces to KG equa- 
tion, provided that the angular momentum quantum number is given through 
the relation / = j ± 1/2, and the radial KG wave function is the radial large- 
component of the Dirac spinor. A parameter A = 0, 1 has been introduced 
in the spin-dependent term U{r) so that when A = the Dirac equation 
reduces to the KG equation. The formalism developed in Sec II has been in- 
vestigated through the 4-vector Coulomb, the scalar linear, and the equally 
mixed scalar and 4-vector potentials in Dirac and KG equations. SLET 
results were found in excellent agreements with the exact and numerically 
predicted results found elsewhere. 

To the best of our knowledge, this is the first time that both Dirac and 
KG equations, with arbitrary scalar and/or 4-vector potentials, have been 
solved through one analytical procedure that yields highly accurate and fast 
converging results. 

In conclusion, SLET is to be understood as being an expansion through 
some existing quantum number, which depends on the symmetry of the po- 
tential of interest, in any Schrodinger - like eigenvalue problem. For example, 
in the case of spherically symmetric potentials the expansion parameter is 
l/l where I = I — j3 and in the case of cylindrically symmetric potentials 
I = \'m\ — j3 where m is the magnetic quantum number. This understanding 
suggests we should not worry about the A^-dimensional form of the wave equa- 
tion in hand and we should expand directly through the quantum numbers 
involved in the problem. Therefore, the difficulties associated with inflat- 
ing the dimensions of the Dirac equation are resolved by using SLET. The 
attendant method seems more flexible than the SLNT. 
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Finally, we would like to mention that the formalism we have introduced 
in this work is also applicable to more complicated relativistic potentials. In 
particular, screened Coulomb potentials which have great utility in atomic, 
nuclear, and plasma physics. Also, the wave functions, and normalizations for 
spin-1/2 and spin-0 particles in spherically symmetric scalar and/or 4- vector 
potentials can be obtained. 

Appendix 

The definition for 5's and 5[s in ai and 0:2 of Ref. [1] are as follows 



£1 = -2(2/3+1) , £2 = 3(2/3 + 1) 
^3 = -4 + ^[7"'(ro) + 2\/"'(ro)i?o" 
<;i = 5 + ^[7""(ro) + 2V'""(r„)B„ 



82 = 3/3(/3 + 1) + 



Q 

rtV"{ro)E2 

Q ' 



<53 = -4(2/3+1) , (54 = 5(2/3+1), 



7 {ro) + 2V {r-o)Eo 



i2og 



5fi = 7 



720Q 



7 (r„) + W {r;)E, 



(31) 
(32) 

(33) 

(34) 

(35) 
(36) 
(37) 

(38) 



The terms including Ex have been dropped from the expressions above 
since -Ei=0. 
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Table 1: Comparison of the energy levels (in keV units) for spin-0 particle 
in V{r) = —a/r and S{r) = 0. W represents the binding energies, where 
Wo = (Eo - l)m^c\ W2 = {Eo + E^/f - l)m^c\ and W^ = {E^ + E^/f + 
E^/l^ - l)m^c\ 



State 


-Wo 


-W2 


-W^ 


Exact [26] 


Is 


3.71649 


3.71653 


3.71656 


3.71659 


2s 


0.92910 


0.92910 


0.92911 


0.92912 


2P 


0.92910 


0.92910 


0.92910 


0.92909 


3s 


0.41293 


0.41293 


0.41293 


0.41294 


3P 


0.41293 


0.41293 


0.41293 


0.41293 


M 


0.41293 


0.41293 


0.41293 


0.41293 


As 


0.23227 


0.23227 


0.23227 


0.23228 


AP 


0.23227 


0.23227 


0.23227 


0.23227 


Ad 


0.23227 


0.23227 


0.23227 


0.23227 


4/ 


0.23227 


0.23227 


0.23227 


0.23227 
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Table 2: Comparison of the energy levels (in eV units) for spin-1/2 particle 
in V{r) = —a/r and S{r) = 0. W represents the binding energies, where 



Es/P 



{Eo - l)mec2, 



Wo 



{Eo + E^/f - l)m,c\ and W3 = {E, + E^/l^ + 



State 


-Wo 


-W2 


-w^ 


Exact [26] 


lsi/2 


13.60639 


13.60657 


13.60667 


13.60603 


2Sl/2 


3.40151 


3.40153 


3.40155 


3.40152 


2P1/2 


3.40150 


3.40151 


3.40151 


3.40152 


2P3/2 


3.40147 


3.40147 


3.40147 


3.40148 


3Si/2 


1.51177 


1.51177 


1.51178 


1.51178 


3Pl/2 


1.51177 


1.51178 


1.51177 


1.51178 


3P3/2 


1.51177 


1.51177 


1.51177 


1.51177 


34/2 


1.51177 


1.51177 


1.51177 


1.51177 


30?5/2 


1.51176 


1.51176 


1.51176 


1.51177 


4Si/2 


0.85037 


0.85037 


0.85037 


0.85038 


4Pl/2 


0.85037 


0.85037 


0.85037 


0.85038 


4P3/2 


0.85037 


0.85037 


0.85037 


0.85037 


44/2 


0.85037 


0.85037 


0.85037 


0.85037 


4(^5/2 


0.85037 


0.85037 


0.85037 


0.85037 


4/5/2 


0.85037 


0.85037 


0.85037 


0.85037 


4/7/2 


0.85037 


0.85037 


0.85037 


0.85037 
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Table 3: KG results for part of the mass spectra (in GeV) of ^, ^ system 
with S{r) = Ar, A = 0A37GeV^ and m = lA2GeV. The values in round 
parentheses are those of SLNT [6]. The values in square parentheses are 
those ofRef. [17]. 



1 


rir 


1 


2 


3 


3.12 


3.46 


3.74 


3.99 


(3.12) 
3.1] 
1 3.7 


(3.46) 

[ 3.47] 

3.95 


(3.74) 
3.73 

4.18 


(3.99) 

3.98 

4.387 


(3.7) 
3.7 
2 4.15 


(3.96) 
3.95 
4.36 


(4.18) 
4.17 
4.56 


(4.387) 

4.387 

4.74 


(4.16) 

4.17 

3 4.537 


(4.36) 
4.38 
4.72 


(4.56) 
4.56' 
4.89 


(4.74) 
4.73 
5.05 


(4.537) 
4.54 
4 4.88 


(4.72) 
4.72 
5.04 


(4.89) 

[4.9] 

5.195 


(5.06) 
5.05 
5.35 


(4.9) 

[4.89] 


(5.04) 

[5.04] 


(5.198) 
5.195 


(5.35) 
5.35 
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Table 4: Dirac results for part of the J/^ mass spectra (in GeV) with S{r) = 
Ar, A = 0.137 GeV^, m=1.12 GeV, and j = 1 + 1/2. The values in round 
parentheses are those of SLNT [12]. The values in square parentheses are 
those ofRef. [17]. 



1 


rir 


1 


2 


3 


3.094 


3.432 


3.712 


3.958 


(3.106) 


(3.449) 


(3.733) 


(3.982) 


3.103 


3.442 


3.725 


3.973 


1 3.691 


3.935 


4.158 


4.363 


(3.698) 


(3.948) 


(4.175) 


(4.383) 


[3.7] 


3.946 


[4.170 


4.377 


2 4.146 


4.348 


4.538 


4.718 


(4.152) 


(4.359) 


(4.553) 


(4.736) 


4.158 


[4.36] 


4.551 


4.732' 


3 4.531 


4.707 


4.876 


5.037 


(4.535) 


(4.717) 


(4.889) 


(5.054) 


'4.545' 


4.72 


4.89 


5.053 


4 4.870 


5.03 


5.183 


5.331 


(4.874) 


(5.038) 


(5.195) 


(5.345) 


'4.886 


5.043 


5.198 


5.346 
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Table 5: Dirac results for part of J/^ mass spectra (in GeV) with S{r) = Ar, 
A = 0A37GeV^ and m = lA2GeV, and j = I — 1/2. The values in round 
parentheses are those of SLNT [12]. The values in square parentheses are 
those ofRef. [17]. 



1 


rir 1 


2 


3 


4 


3.4726 


3.7611 


4.0118 


4.2369 


(3.471) 


(3.760) 


(4.010) 


(4.236) 


3.47 


3.757 


4.006 


4.23 


1 3.9646 


4.1961 


4.4069 


4.6020 


(3.954) 


(4.194) 


(4.406) 


(4.600) 


3.965 


4.194 


4.403 


4.597 


2 4.3713 


4.5699 


4.7552 


4.9295 


(4.370) 


(4.568) 


(4.754) 


(4.920) 


■ 4.374] 


4.560 


4.753 


4.926 


3 4.7263 


4.9030 


5.0700 


5.2290 


(4.726) 


(4.902) 


(5.069) 


(5.228) 


4.731 


4.95 


5.07 




4 5.0456 


5.2062 


5.3595 


5.5065 


(5.045) 


(5.205) 


(5.359) 


(5.507) 


5.053 


5.21 


5.361 


5.506 



Table 6: Spin- averaged Dirac eigenvalues e„^i, Eq(30), for the equally mixed 
scalar and 4-vector power-law potential V{r) = S{r) = Ar^ + Vo-, Vo = 
-2.028GeV, v = 0.1, A = a''+\ a = 1.709, and m^ = 1.6179^6^. 



cc 


rir 


1 


Ref.[19 


SLNT [5] 


SLET 








1.2364 


1.2400 


1.2358 


1 





1.3347 


1.3400 


1.3343 


2 





1.3923 


1.3980 


1.3915 


3 





1.4335 


1.4390 


1.4326 


4 





1.4657 


1.4710 


1.4648 





1 


1.3071 


1.3090 


1.3072 


1 


1 


1.3731 


1.4110 


1.3730 





2 


1.3544 


1.3580 


1.3540 



Table 7: Spin-averaged Dirac mass spectra (in Gev units), M„^/ = 2E, for 
the bb system in the equally mixed scalar and 4-vector power-law potential 
of Eq.(29). 



rir 


1 


Ref.[19 


SLET 








9.4390 


9.4347 


1 





10.0260 


10.0228 


2 





10.3694 


10.3640 


3 





10.6147 


10.6084 


4 





10.8063 


10.800 





1 


9.8613 


9.8612 


1 


1 


10.2549 


10.2537 





2 


10.1435 


10.1406 
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Table 8: Spin-averaged Dirac mass spectra (in Gev units), M„^; = 2E, for 
the ss system in the equally mixed scalar and 4-vector power-law potential 
of Eq.(29). 



fir 


1 


Ref.[19 


SLET 








1.0040 


0.9999 


1 





1.5532 


1.5500 


2 





1.8804 


1.8752 


3 





2.1162 


2.1099 


4 





2.3012 


2.2950 





1 


1.3976 


1.3974 


1 


1 


1.7710 


1.7697 





2 


1.6647 


1.6619 



Table 9: Spin-averaged Dirac mass spectra (in Gev units), Mn^i = Ed + rric, 
for the cd system in the equally mixed scalar and 4-vector power-law potential 
of Eq.(29). 



Ur 


1 


Ref 19 


SLET 








1.9730 


1.9713 


1 





2.2107 


2.2093 


2 





2.3601 


2.3577 


3 





2.4699 


2.4671 


4 





2.5571 


2.5542 





1 


2.1413 


2.1413 


1 


1 


2.3097 


2.3092 





2 


2.2612 


2.2599 
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